On the formula of regularized traces 



o 



> 



O 



A.I. Nazarov, D.M. Stolyarov, P.B. Zatitskiy* 
January 18, 2013 

We consider an operator L generated by a 2m-order differential expression 



• 2m-2 

CM : ly^ {-iTD'^^y + ^ Vk{x)D''y (1) 

I— !■ 

C/^ \ and by m boundary conditions 

'S; P,(D)y(0) = 0, j = l,...,m. (2) 



Here pk ^ -^i,«oc(IR+) are real functions while Pj is a polynomial of degree kj\ moreover, 
the system of boundary conditions is assumed normalized, i.e. Q < ki < k2 < ... < < 
2m - 1. 



]|0 ■ Suppose this operator is self-adjoint in the space L2(M_|_), semibounded from below 

and has purely discrete spectrum {A„}^^ enumerated in ascending order according to 
■ the multiplicity. 

! Let Q be an operator of multiplication by a real function q G Loo(IR+)- Then an 

operator L + Q also has a purely discrete spectrum {fin}t=i- 

A number of papers beginning from the pioneering work are devoted to the cal- 
culation of spectral functions and regularized traces of differential operators. The aim of 
^ ! our paper is to prove the following statement. 



Theorem 1. Let q have a bounded support, and the function tjj{x) = - J q{t)dt has 



hounded variation at zero. Then the following relation holds: 



oo 



5, ^ - A„ - ^ / qm] = -H0+) ■ (y - i - ^ (3) 

n=i 



where 







1 1 



ci = Xl"" ; Cn = A^™ - A^™i, n > 1; x = XI 
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^This condition is used only to conclude the relation from Proposition [TJ It seems that this 
condition may be weakened. To do this one needs to obtain a global estimate of the difference of spectral 
functions to operators L and L. However, up to the moment we do not know such an estimate. 



Remark. Formula (|3]) was conjectured by A.I. Nazarov at the conference in Moscow, 
2007, during the A.S. Pechentsov lecture. For one-term boundary conditions Pj{D) = D^^ 
this formula was proved in the preprint pjl. 

We introduce two auxiliary operators: 

L is a 2m-order operator with boundary conditions ([2]), self-adjoint in L2(M+); its 
lower-order coefficients are compactly supported and coincide with lower-order coefficients 
in ([1]) on the segment [0, i?] D supp(g); 

Lq is the operator (— l)"^/}^'" with boundary conditions (IJjJfl. 

For A G M we denote by 6'(x, A) the spectral function of the operator L, i.e. the kernel 
of its spectral projector Ex, see [6j. In a similar way, 9{x,y, A) is the spectral function of 
the operator L. Also we denote by Ho{x, y, r) the Green function of the operator Lq — r 
(for r ^M+). 

Let C be the value of such that arg(^) e [0, — [, while z = exp(— ). We introduce 
the matrix B{C) = [Pjiiz^'K)]^^^^ and define A(C) = det(i3(C)). 

The following statements were proved in see also [S]. We state them with some 
redefinitions. 

Proposition 1. A (^7, Ch. 1, Sec. 2]). Let the following condition hold: 

A. The matrix B{() is non- degenerate, and the elements of the inverse matrix satisfy 
the estimate [-6-^0]^^ = 0(|Cr^O as |CI ^ oo- 

Then the relation 



holds uniformly on [0,-R]^. 

First, we note that he assumptions of Proposition [T] may be weakened. 

Lemma 1. For \(\ large enough, condition A is always satisfied. Thus, this condition in 
Proposition [I] can be omitted. 

Proof. We claim that A(C) is a x-degree polynomial of (. Indeed, the columns of B{Q are 
fcj-degree polynomials, while the highest degrees cannot reduce since the corresponding 
coefficient is the Vandermonde determinant det{W{z''^ , . . . , z^"^)) ^ 0. Thus, A(^) ^ 
for sufficiently large 

Further, the cofactor of is evidently a polynomial of C, and its degree does not 
exceed k —ky The statement now follows from the Cramer formula. □ 

^ Three particular cases: 1) kj = 2j — 2; 2) kj = 2j — 1; 3) kj = j — 1 were considered earlier in papers 




as A — > +00, 



holds uniformly on R_|_. 

B (j^ Ch. 1, Sec. 4])- Let condition A be satisfied. Then the relation 



9{x,y,X) = 9{x,y,\) + o(l). 



as A — )■ +00 




Proof of Theorem Ql Note that the left-hand side of ([3]) can be rewritten as follows, see O 
proof of Theorem 1]: 



hoo 



q{x) ■ I 9{x, X, A) ) dx. 



By Lemma [1] this formula can be transformed to 

Si = lim / q{x) ■ ( / HQ{x,x,T)dT ] dx. (4) 

A-i.+oo J \ 2ni J IT / 

|t|=A 

Now we use the explicit formula for Hq obtained in P| Lemma 1] . We state them with 
some redefinitions. 

where the determinant Aap{() is obtained from A(^) if we substitute Pj{—iz°'~^() for 
Pj{iz^~^C) in /3-th line. Surely, for this formula we need A((^) 7^ but this is the case for 
|r| large enough. 

We change the variable r = (^^™ in the inner integral in (j4j). It is easy to see that 

idr f d( 



r tdr v^^.-i^_ f 
J 2mC2™-i ^ J 



2711 J 2mC2™"^ ^ J -nil- z) vr 

|r|=A "-^ Fa 

(here Fa is the arc of the circle {C = Aa^e*"^ : g]0, ^[}). Therefore (jl]) is rewritten as 
foUowfl: 

+00 



r 



A 



We denote B„« = lim ^f^/P and claim that 



-\-oo 

Si = ^ lim / q{x) ■ [ V exp(^(z"-i + z'^-^) <) ■ dC dx. (6) 



o,/3=l 

Indeed, since Aq,^ is a polynomial of ( and its degree does not exceed x, the relation 

U + OiC) (7) 



A(C) 



^In [3J Theorem 2] a formula similar to (O was obtained under additional assumption: A(C) ^ for 
r £ R \ {0}. This formula contains the real part of the integral over the segment [0, A^^] instead of 
integral over T\. In fact, this formula is correct only under a stronger assumption: =^ for r ^ 0. 

In this case [H formula (0.5)] is reduced to (O by the Cauchy theorem. 



holds as ICI — )■ oo. Further, for any 1 < a, P < m the inequahty < aTg{z°' ^ + ^) < 
^^^TT holds true. Therefore, the integrals 



are bounded by ([7]) uniformly with respect to x > and tend to zero as x > 0, A — )■ oo, by 
the Jordan Lemma. By the Lebesgue Dominated Convergence Theorem, lim J q{x) ■ 

A— >+oo Q 

Xq,^(x) dx = 0, and (El) follows. 

We calculate the inner integral in (E]) and rewrite the formula for Si as follows: 

oo 

Si = ■ lim / q{x)g{sx)s dx, (8) 

2711 s^oo J 


where 

f\ TO m exp{i{z" + z^) y) - exp{i{z''-^ + z^^-^) y) 

9\y)= 2^ ^I3a^ap , 

a,l3=l ^ 

while 

1 



To pass to the limit in ([8]) in terms with a = P = 1 and a = P = m, one uses the 
Riemann locahzation principle, see, e.g., Ch. I, Sec. 33], and the Vallee-Poussin test 
of convergence, see [TUl Ch. Ill, Sec. 3]. Other terms, after integration by parts, are 
covered by the Lebesgue Theorem, since the exponents have negative real parts. Thus we 



arrive a^^' 



oo 



5i = ^ ■ ^(0+) I g{y) dy = ■ Sp(PB) (9) 



(the last equality follows from [HI 3.434.2]). 

Note that the entries of matrix B are quotients of determinants composed of the 
leading coefficients of polynomials that are entries of determinants Aa/siC) and A{(). 
Direct calculation via the Cramer formula gives 

B = W ■ diag[(-l)'=\ . . . , (-1)'=™] • W"^ (10) 

(we recall that W = (W^j) is the Vandermond matrix generated by the numbers wj = 
j = l,...,m). 

Formula (ITOl) shows that if e^^^ = [l,Wj, . . . , w^''^]'^ is j-th column of matrix W then 
]gg(i) = (^_x)fcjgO)_ Thus, vectors e^^^ form the eigen basis of B. 

Next, by the geometrical progression sum formula, we rewrite the matrix P as follows: 

oo 

P= lim J](-1)V^„C, (11) 
p—>-i— ' 

n=0 

where ^„ = (1, z", . . . , .., ^("^-i)")^. 

Denote by IK the set {kj + 2mp \ 1 < j < m, p E 



'Note that passage from (|8]) to (H)) follows also from [5] Lemma 2] 



Lemma 2. If n e K then Sp((^„^^B) = (-l)'^m. If n ^ K then Sp((^„^^B) = 
(-l)"+im. 

To prove this lemma we need an obvious proposition. 

Proposition 2. Let u, v be elements of Hilbert space and let a,b E C. Then {au + bv, u) = 
implies {au + bv, —au + bv) = \au + top. 

Proof of Lemma\^ First, 

We consider two subspaces 

U = Lin{e^^^ \ kj = 1 (mod 2)}; V = Lin{e^^^ \ kj = (mod 2)}. 

Note that U + V = C™, since the vectors e^^^ form a basis. Thus, there exists a 
decomposition (fn = u + v with u G U, v G V. Further, we have shown that B|^ is an 
identity operator while B|^ multiplies by —1. This implies M(pn = —u + v. 

1. Let n G K. Then is an eigenvector of B corresponding to the eigenvalue (—1)"^. 
Therefore (B(^„,(^„) = (-1)"((^„, = (-l)"m. 

2. Let n ^ IK and n = 1 (mod 2). Then direct calculation shows that ± f/. By 
Proposition m (M(pn,fn) = \ fn\^ = fn- 

3. In a similar way, if n ^ IK and n = (mod 2) then (f V. By Proposition [21 
= |v?„P = m. □ 



To complete the proof of Theorem we define the set K = {ki, k^}- Then formula 
TT]) and Lemma [2] imply 



n=0 



Sp(PB) = hm ^(-l)>"Sp(^.CB) = 

oo oo 
jeK p=0 0<j<2m,j<^K p=0 



777- 

ie-fsT 0<j<'2m,jfK 



Since the set i^' contains exactly m elements, we get 

-I (>> - I 



Sp(PB)=rn^lim (J^i^- 5: 



jeK ~ ^ 0<j<2mJ^K ^ 



^lE^- E ^) = ^ E ^-E^' 

ieX 0<j<2m,jfK 0<j<2m jG-ftT 



m{2m — 1) , , , 
(12) 



Theorem follows immediately from (|9]) and (fT2|) . □ 
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